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We study the entanglement of families of Unruh modes in the Bell states I'l’*) = l/\/2(|00) ± 111)) 
and 14/^) = l/%/2(|01) ± |10)) shared by two accelerated observers and find fnndamental differences 
in the robustness of entanglement against acceleration for these states. States are entangled for 
all finite accelerations, whereas, due to the Unruh effect, states lose their entanglement for finite 
accelerations. This is true for Bell states of two bosonic modes, as well as for Bell states of a bosonic 
and a fermionic mode. Furthermore, there are also differences in the degradation of entanglement for 
Bell states of fermionic modes. We reveal the origin of these distinct characteristics of entanglement 
degradation and discuss the role that is played by particle statistics. Our studies suggest that the 
behavior of entanglement in accelerated frames strongly depends on the occupation patterns of the 
constituent states, whose superposition constitutes the entangled state, where especially states 4?^ 
and exhibit distinct characteristics regarding entanglement degradation. Finally, we point out 
possible implications of hovering over a black hole for these states. 
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I. INTRODUCTION 

Entanglement and quantum correlations in general 
play an important role in different areas of physics, as, 
for example, in quantum information [l[ and black holes 
ii- Furthermore, it is known that the quantum cor¬ 
relations of an entangled state shared by accelerated ob¬ 
servers are not invariant with respect to acceleration but 
are altered by the Unruh effect 0. Interestingly, the 
noninvariance of quantum correlations in this relativistic 
regime can be employed to carry out quantum informa¬ 
tion tasks il. Although accelerated motion can, in 
some special cases, create entanglement between Unruh 
modes 0, generally entanglement is degraded due to 
the Unruh effect. 

In the past, the degradation of entanglement in bi¬ 
partite states composed of Unruh modes shared by an 
inertial observer and an uniformly accelerated one was 
studied in detail [TTI - [^ . For fermionic fields, entangle¬ 
ment approaches a finite value in the infinite acceleration 
limit , while for bosons it vanishes asymptotically Q . 
One way to study the entanglement between two accel¬ 
erated observers is to analyze an entangled state shared 
by three parties, where two parties are in accelerated 
motion, and subsequently trace out the inertial observer 
[234^ . A more natural way to study entanglement in 
this framework is to restrict to bosonic entanglement in 
accelerated two-mode squeezed states and use tools from 
continuous variable quantum mechanics . There 

it was found that, for these squeezed states, entangle¬ 
ment vanishes for finite accelerations, in contrast to the 
entanglement of states shared by an accelerated observer 
and an inertial one. 

As realized more recently in [29l - l3lj| . there are some 
caveats in the interpretation of states of Unruh modes. 


Still, however, the use of Unruh modes, which allows for 
closed analytical solutions, provides a valuable frame¬ 
work to understand the mechanisms that lead to a de¬ 
crease of quantum correlations in entangled states when 
described by accelerated observers. The goal of this work 
is to provide further insight into the degradation of en¬ 
tanglement that occurs when entangled states are ob¬ 
served by accelerated parties. Therefore, in this work, 
we study families of states composed of Unruh modes 
that are maximally entangled from the inertial perspec¬ 
tive and investigate the residual entanglement when these 
states are seen by uniformly accelerated observers. 

We start by studying the entanglement between 
two accelerated observers sharing the fermionic Bell 
states !<!)='=) = l/A/2(|0cj0n) ± |lcvln)) and = 

l/-\/2(|0tjln) ± |lw0n)). We find that entanglement is 
nonvanishing for all accelerations, and that the degra¬ 
dation of quantum correlations depends on the specific 
state shared by the parties. The reason for the survival 
of entanglement for fermions is rooted in the statistics 
obeyed by fermions, while the statedependence of entan¬ 
glement degradation seems to be originated in the occu¬ 
pation pattern of the state that is shared between the two 
observers. By occupation pattern we mean the pattern of 
both constituent states (in the following just called con¬ 
stituents), for example, |00), |01), |10), and |11), whose 
superposition defines the entangled state. That is, the 
set of excitations created by the Unruh effect depends on 
the state one starts with. 

Furthermore, we study maximally entangled bosonic 
states shared by two accelerated observers and, in con¬ 
trast to [^, , we consider two bosonic modes of 

frequency uj and U, respectively, in the Bell states 
and We find that the Unruh effect degrades entan¬ 
glement in these two states very differently. While, as 


2 


in [13, m, 113, bosonic modes in state lose all their 
entanglement for finite accelerations, the entanglement 
of these modes in state is nonvanishing for all finite 
accelerations. We find that this crucial difference in the 
degradation of entanglement is due to the differing occu¬ 
pation patterns of the constituents of the two Bell states, 
and is manifest in the appearance of a “cut-off function” 
in the expression of the negativity. 

Then we extend our studies to accelerated states of a 
bosonic mode maximally entangled with a fermionic one. 
Thereby we find states whose negativity factorizes. How¬ 
ever, due to the bosonic mode involved in these states, 
there is no entanglement surviving in the infinite accel¬ 
eration limit. Then, moving on to Bell states and 
we obtain qualitatively the same behavior as for the 
purely bosonic Bell states. Thus, we find evidence for 
the importance of the occupation patterns of the con¬ 
stituents. The particular occupation patterns of its con¬ 
stituents can protect a state’s entanglement against the 
effects of acceleration, i.e., the Unruh effect. Using an ef¬ 
fective state picture we are able to explain the differences 
in the behavior. 

In the past, the different behavior of the entanglement 
of accelerated fermions and bosons led to some discus¬ 
sions Here we address this issue and discuss 

the role played by particle statistics by combining the re¬ 
sults of the fermion-fermion, boson-boson and the boson- 
fermion cases. We conclude that there are essentially two 
factors determining the fading of entanglement, where 
one is purely from particle statistics and one strongly de¬ 
pends on the occupation patterns of the constituents of 
the state that is considered. Finally, we discuss some ef¬ 
fects that hovering over a black hole at a fixed distance 
from the horizon has on the entanglement of the states 
that are studied in this work. 

The outline is the following. In Sec. [Ill we give a short 
introduction to quantum fields in Rindler space. In Secs, 
cniiia and El we study the entanglement of fermion- 
fermion, boson-boson and boson-fermion Bell states in 
accelerated motion, respectively. Then, in Sec. ED we 
discuss the role of particle statistics for entanglement 
degradation and give two factors that determine the char¬ 
acteristics of said degradation. In Sec. IVHI we outline 
possible implications of our findings for Bell states in the 
vicinity of a black hole and, finally, in Sec. IVHII we give 
the conclusions of this work. 

For the sake of brevity, throughout this work, we call 
the occupation patterns of the constituents of a state just 
the structure of a state. 


II. QUANTUM FIELDS IN RINDLER SPACE 

We give a brief introduction to quantum field theory 
in Rindler space. The main purpose is to introduce the 
framework we are using in this work. More details can 
be found, for example, in [3, (H, We work in units 
where c = h = ks = I- Minkowski coordinates {t,x) and 



FIG. 1: (Color online) Rindler space: Regions / and II 
are causally disconnected due to the presence of horizons 
at X = t and x = —t. The worldlines of two observers A 
and B in (differently) accelerated motion are shown. The 
states shared by A and B are prepared and distributed 
in the asymptotic past. 


Rindler coordinates (^, rf) are related by the transforma¬ 
tions 

t=^sinh( 77 ), (la) 

a; =^cosh(? 7 ), (lb) 

where the range of ^ and r] is given by —oo < ^,7] < oo. 
Notice that ^ is positive in the right wedge (region I) 
and ^ is negative in the left wedge (region II). Then we 
obtain the following metric 

ds^ = dr^ — . ( 2 ) 

Considering a world line with ^(r) = i, where r is the 
proper time along this trajectory and |a| is the proper 
acceleration, we find rjij) = ar. Thus, in Minkowski co¬ 
ordinates the world line reads ^(t) = isinh(ar), x{t) = 
icosh(aT). As shown in Fig. (T] the two regions I and 
// of Rindler space are causally disconnected due to the 
presence of horizons at x = t and x = —t. A timelike 
Killing vector in region / is given by 9^ i~dri in II). 

A. Bosons 

We consider the quantization of a massless scalar field 
(j) (see m, 113 for details). We quantize fields with re¬ 
spect to the Killing vectors 9^ and —dri independently in 
the two regions. The Klein-Gordon equation □(/> = 0 in 
Rindler coordinates has solutions which depend on 77 as 
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where a) is a positive parameter and the sign ± depends 
on the Rindler wedge. These modes are positive fre¬ 
quency modes with respect to the respective tiinelike 
Killing vectors. We denote by i.e., oc 

and ujjj oc the solutions in regions I and II, re¬ 

spectively. We can expand (j) in this basis and obtain 


(j) = 




aiui + a^ui* 




*), (4) 


where the and are the usual commuting an¬ 

nihilation and creation operators in regions I and II, 
respectively. The dependence on the proper time t = ^ 
is given by oc Therefore, the energy uj seen by 

the accelerated observer is given by a; = aw. Remember 
that in Minkowski space (j) can be expanded as 


poo 

(j)= duJM + 

Jo 


Mt M*\ 


(5) 


where and are the commuting Minkowski an¬ 
nihilation and creation operators. These two expansions 
lead to different Fock spaces. Consider the Minkowski 
(M) vacuum |0 )m and the Rindler (i?) vacuum |0)fi = 
| 0 )/ 0 | 0 )// that are defined as 

a^| 0 )M= 0 , ( 6 a) 

ai\0)R=ai^\0)R = 0, ( 6 b) 


and, in general, a^’'^'^|0)M 7 ^ 0. Next, we want to intro¬ 
duce the so-called Unruh basis that we use in the fol¬ 
lowing. The Rindler creation and annihilation operators 
are related to the corresponding Unruh ([/) ones 
i^y Bogoliubov transformation as follows [3^ : 


a,', = 


alf = 


1 


sinh(7rw) 

1 

y/2 sinh(7rw) 


e 2 a, 


(7,2 , - 21 ^ (7,lt\ 

'£D +e ^ a- 'j. 

(7a) 

(7,1 , -2L^ (7,2t\ 

+e ^ a-’ '1 . 

(7b) 


These share the positive frequency analyticity properties 
of the u/'f and therefore have the same vacuum state 

UJM 

4 ’'’>W = 0 , i.e., |0 )c/ = | 0 )m M- 
In this work we study accelerated observers confined 
to Rindler wedge I and work in the Unruh basis. To ob¬ 
tain the appropriate description of what an accelerated 
observer is experiencing, we use Bogoliubov transforma¬ 
tions © to go from the Unruh basis to the Rindler basis. 
For an Unruh mode w, the vacuum and one particle states 
are given by 


I \ \^ tanh (r*). . . . 

\0c.)u =2^-—^\na,)i\nQ)ii, 

ii \ tanh (r) , -—. -, \ i \ 

\^Cj)u = / - .2/ X V ^ + l|n + 

cosh (r) 

n ^ J 


( 8 a) 

( 8 b) 


for a massless uncharged scalar field, where the acceler¬ 
ation parameter r is set by w and they are related by 
r = arctanh(e“'^‘^). Further, |n)/ and \n)ii are the n- 
particle states in regions I and II, respectively. Note 
that in (iRbl) we denote the one-particle state a^’^^ | 0 (i)c/ 
by |lii)f 7 . That is, the excitation is localized in I. Simi¬ 
larly, for a massless charged scalar field the vacuum and 
one-particle states are given by [l^ 


, , v-tanh"+'"(r), , , , 

|0 £d)(7 =2^-|nTO)j|nTO)j/, 


cosh (r) 


(9a) 


, tanh”+'"(r) , -- .n \ i \ 

^)ij = V- , 3, ' yn -h l|(n -h l)m)7|nm)77, 

cosh (r) 

n.m ' I 


(9b) 


^tanh"+™(r) , -, ,,, , , 

| 1 cd ) c / = > - --^--—\/m + l\n[m + l))i\nm)u, 

cosh (r- 

n,.m. k I 


(9c) 


where \nm)ii denotes the state of m/n particles/ anti¬ 
particles of energy w in region R = I, II. 


B. Fermions 


We model the fermionic field by a massless Grassman- 
nian valued scalar field ijj. Then the quantization of ^|J 
can be carried out analogously to the bosonic case. To 
obtain the appropriate description of what an accelerated 
observer is experiencing, we use the Bogoliubov transfor¬ 
mations 0 
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d! = 


■\/2cosh(7rw) 

1 

a/ 2 cosh(7rw) 


( 21 ^ u ,2 , - 

e 2 a.r + e 






..) 


( 10 a) 

( 10 b) 


to go from the Unruh basis to the Rindler basis. We 
choose the notation \ijkl)a, = \ia,)^ 0 |jii )77 G \ka,)~[ 0 
where +/— denote particles and anti-particles, re¬ 
spectively. We use the chosen order throughout the work 
and obtain for an Unruh mode w 


|0i )(7 = cos^(r/)|0000)^ - cos(r/) sin(r/)|0011)(i, 

-|-cos(r/)sin(r/)|1100)£D - sin^(r/)|llll)£D, 

(11a) 

\^q)u =cos(r/)|1000)£D - sin(r/)|1011)£D, (11b) 

\^E)u =cos(r/)|0010)£D -f sin(r/)|1110)£D, (11c) 

where the acceleration parameter rj is given by ry = 
arctan(e“’^‘^). 


C. Unruh effect 

An accelerated observer does not necessarily agree with 
an inertial observer on the number of particles in a given 
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state. Consider, for example, the Minkowski vacuum 
state |0)m = |0)(7 = |0). Then the vacuum expectation 
value of the (Rindler) number operator 
can be calculated using dini), leading to 

(0|ar'an0) = l + e-. (12) 


Thus, we see that an accelerated observer perceives the 
Minkowski vacuum as a thermal state of temperature Tjj 
(Unruh temperature) [ssj. 


Tu = 


a 


(13) 


Further, we introduce the fermionic and the bosonic par¬ 
tition functions Zp and Zp that are given by 


ryijJ _ 

— 


Z^p=\ 


1 — e 
e~ 


1 


1 — e 
= 1 -I- e“ 


(14) 

(15) 


There are some subtleties when working with the 
global Unruh modes (0), di), and (HD, as pointed out 
in [ 29 I - I 31 II . Therefore, one point that we want to empha¬ 
size is the implicit dependence on the acceleration a in ([7|) 
and m through the relation uj = —. As a consequence, 
after fixing the frequency w, each of the Unruh modes 
(©,(191), and (ED forms a family of modes that is labeled 
by a. That is, by varying the acceleration parameter r/rj 
(or equivalently a), one also varies the particular Unruh 
mode under consideration. In order to pick a particular 
state, Lo and a have to be fixed. Intuitively, one can say 
that the acceleration a is already encoded in the Unruh 
modes. We revisit these issues when we discuss the en¬ 
tanglement of fermions in Sec. Hill In the following, for 
the sake of simplicity, we omit the tilde in w whenever it 
is clear from the context whether we are talking about w 
or UJ. 

In this brief introduction to quantum fields in Rindler 
space we have set up the tools and notation we are us¬ 
ing in the following. In the next section we study the 
degradation of entanglement in fermionic Bell states due 
to acceleration. 


III. ENTANGLEMENT AND ENTROPY OF 
UNIFORMLY ACCELERATED FERMION 
STATES 


Due to the anticommutativity of fermionic creation 
and annihilation operators (Pauli exclusion principle), 
there is only a finite number of maximally entangled 
states of two modes of a fermionic field. Considering par¬ 
ticle states, there are just two possible maximally entan¬ 
gled states. These are the two Bell states of two fermionic 
modes (FF): 

=± {\l!;)^\0E)u ± |05)c/|l^)^) , (16a) 

\<^$p) =± {\o!;)uK)u ± lO^lla)^) ■ (16b) 


The subscript U emphasizes that we are working in the 
Unruh basis. 


A. Negativity 


We consider two families of entangled fermionic Un¬ 
ruh modes uj and D undergoing constant accelerations 
Guj and an, respectively. The acceleration parameters of 
the modes are denoted by rj and . Therefore, starting 

from the families of states {tpi} = written 

in the Unruh basis, we use ED to obtain the density 
matrices Pj// 

P^iji = ( 17 ) 

To describe the system as it is seen by an observer con¬ 
fined to region /, we have to trace out modes that have 
their support in the inaccessible region. Then the re¬ 
duced density matrix pi is given by 

p, = Trn [pfj!) ■ (18) 

As a measure of entanglement we use the negativity N 
(defined in Appendix]©, which is an entanglement mono¬ 
tone [ 3 ^. It is known that a bipartite state is not sep¬ 
arable if its negativity is nonzero [s^. Furthermore, the 
vanishing of the negativity provides a necessary and suf¬ 
ficient condition for the separability of mixed states of 
two qubits [ 3 ^. Fermions, in general, cannot be treated 
as qubits. However, when charge superselection is re- 
spected, two fermionic modes can be represented as two 
qubits [39j. A further property of the negativity is that 
a state with vanishing negativity contains no distillable 
(free) entanglement, although, in this case, there can be 
nondistillable (bound) entanglement present . In this 
work, we ignore the possibility of bound entanglement 
and refer to free entanglement as entanglement. 

To obtain the negativities W of states tpi^ we have to 
calculate the partially transposed reduced density matri¬ 
ces We find that these matrices are block diagonal. 
More details of the calculations can be found in Appendix 
E The final results for the negativities Ni read 



^FF 


1 1 1 
Z Zj p Zj p 


_ rtf + ng j 

(19) 

( 20 ) 


where are the Unruh temperatures (fT51) correspond¬ 
ing to the respective accelerations a^j and an, Zp^^ is the 
partition function (I15II . and w, D are the energies of the 
modes. Further, we introduced the occupation numbers 
and = (1 -|- 

Having obtained the analytic expressions for the neg¬ 
ativities of the families of maximally entangled fermion 
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states dig, we want to comment on the physical inter¬ 
pretation of these states. As discussed in Sec. [Ill we 
are not describing a fixed state ipi, but rather describe a 
two-parameter family of states ipi labeled by and an- 
Therefore, the negativities (IT^ and (l^ give the entan¬ 
glement of the states ipi (maximally entangled from the 
inertial perspective), when the two modes uj and are 
seen by accelerated observers undergoing the accelera¬ 
tions a^, and oq, respectively. Equivalently, we can think 
of states tpi as families labeled by uj and O, when we are 
fixing Uu, and qq. It should be noted that for a given set 
(w, Oij, ao) the only difference between these states is 
the difference in their occupation pattern, in the sense of 
|00) -I- |11) vs. |10) -I- |01). In the following, we discuss 
the effects of acceleration on these families of states and, 
for the sake of brevity, refer to them just as states. 

Considering state it is interesting to note that 

the negativity (1^ factorizes as 

N^±^{rJ,rf) = 2Nf{r-f)Nf{rf), ( 21 ) 


where we denoted = 0) by Nf{rj). Note 

that Nf{rJ) is the negativity in case of only one mode 
being seen by an accelerated observer (acceleration pa¬ 
rameter rj). The negativity Nfir^) = ^cos^ir^) was 
obtained, for example, in M- This product structure is 
absent for state where the negativity is given by 

(O. Thus, the degradation of entanglement in the case 
of a fermionic field shows no universal behavior. Indeed 
different classes of states, and are not equally 

robust against acceleration; see Fig. O That is a feature 
that was absent in previous studies of one accelerated 
observer. 

There is a fundamental difference between states 
and While each state 'ipi is a superposition of 

two states (constituents) of the form \k^)u\lE)u {kj S 
{0,1+}), only for state 'k^^ both such states lead to 
a contribution to a (the same) diagonal element of the 
reduced density matrix that is relevant for entangle¬ 
ment. More precisely, for non-vanishing acceleration, 
| 1 +))}| 0 q )(7 as well as |0+)£/|lQ)y contribute to the 
matrix element |1^)/’|1 q)^(I^|/(IqI/ of the reduced 
density matrix. This is well reflected in the expres¬ 
sion for the negativity, given by (fl^ . The contribu¬ 
tion of |l^)y|0f{)j/ is quantified by np and the one of 
|0+)c/|lQ)y by Due to the symmetry, (fl^ depends 
only the average occupation number rip = l/2{np + n^). 
This behavior distinguishes from 

As in the setting of one mode seen by an accelerated 
observer and one mode seen by an inertial observer fl^ . 

in the limit of infinite acceleration the negativity 
does not vanish, but it approaches a finite limit. The 
surviving entanglement is calculated to be 


lim 

ry ,r)?—J-oo 

lim 

ry ,r)?—>oo 


^ FF 
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(22a) 

(22b) 



FIG. 2: (Color online) Negativities for maximally entan¬ 
gled fermion states {tpi} versus the acceleration a = a^j = 
an, measured in units of (for some length scale L), for 
frequencies oj = ft = jj. For each fixed acceleration a, 
the entanglement degradation for state (blue con¬ 

tinuous line) is stronger than for state 'l'$p (red dashed 
line). The finite asymptotic values for states and 
'^$p are 1/8 and l/4(-\/2— 1), respectively. 


The fact that there is entanglement surviving in this limit 
is specific for initially pure maximally entangled states 
and contrasts with the case of starting from a tripartite 
state, where one observer is inertial and two observers are 
accelerated. In that case, after tracing out the inertial 
observer, the bipartite entanglement between the modes 
observed by accelerated observers vanishes in this limit 
[^ . Numerical studies of fermionic mixed state entan¬ 
glement also showed that entanglement is extinguished 
for most states in the infinite acceleration limit . 

Our results reduce to the known results for one accel¬ 
erated observer if one takes the limit —>■ 0 and we 

obtain the universal behavior reported in [l^ . 

lim Ni = Nf = - cos^(r‘/) (23) 

and thus 

lim Nf = —. (24) 

ry—ioo 4 

Interestingly, the behavior of the negativity under ac¬ 
celeration does not depend on whether there is entan¬ 
glement created in some sectors. We define a sector of 
a state ipi as follows: A sector of state ipi consists of 
all the elements of the reduced density matrix pi that 
contribute to one block of the block diagonal partially 
transposed reduced density matrix For example, 

has four sectors, as can be seen from the partially 
transposed reduced density matrix (see Appendix 

El). When the acceleration is increasing from zero, en¬ 
tanglement decreases in the sector where it is initiated 
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and, depending on the particular structure of the state, 
entanglement is created in previously nonentangled sec¬ 
tors. More details can be found in Appendix lA 21 

The consequences of the fact that states m depend 
on the acceleration only via the ratios Lojuaj and 

n/ao are manifest in (fT^ and (|20ll . One observes that 
high-frequency modes are less effected by acceleration 
than low-frequency modes are. This is due to the larger 
wavelength of low-energy modes. The larger the wave¬ 
length, and therefore the spatial extension, compared to 
the inverse Unruh temperature, the more the system gets 
“stretched” by the acceleration. Thus, the effects of ac¬ 
celeration are stronger in this case and the rate of entan¬ 
glement degradation is higher. 

After this detailed study of the negativity of states 
(in, we now move on to analyze the entropy and the 
mutual information of these states. This provides fur¬ 
ther insight into the effects acceleration has on the cor¬ 
relations in fermion states. 


B. Entropy and mutnal information 

In the following, we analyze entropy and mutual infor¬ 
mation of the fermion states. Since we are considering 
accelerated observers and therefore trace out region II to 
obtain the reduced density matrices pi, the resulting state 
is not pure any more. Thus, the entropy of pi increases 
due to entanglement with modes in region II. Among 
the different measures of entropy, the most widely used 
is the von Neumann entropy S, given by 

S{pi) =-Tri {pi\Yi{pi)). (25) 

In our setting, the von Neumann entropy can be cal¬ 
culated analytically. Since the corresponding expressions 
are quite long and not very enlightening, we give the plots 
of S' as a function of the acceleration in Fig. [3l In the 
limit of infinite acceleration the von Neumann entropies 
approach the asymptotic values S°°(pi), which are given 
in the following 


) = 


ln(32) 3-f2y2, / 3-t-2^2 

In ' 


32 


2^2-3, /3-2V2 

---In 


32 


S-(p.±J= In (8). 


(26a) 

(26b) 


Furthermore, we calculate the mutual information I 
between modes oj and fl as a measure of quantum and 
classical correlations as 


Ii = S {Tr^ {p,)) + S {Tra {p,)) - S{pi), (27) 

where Tr^j^ denotes the trace over mode w/fl. The 
resulting mutual information of states {tpi} (in bits) is 



FIG. 3: (Color online) Mutual information (dashed 
lines), measured in bits, and von Neumann entropy (con¬ 
tinuous lines) for maximally entangled fermion states 
{i^i} plotted versus the acceleration a = a^i = an, mea¬ 
sured in units of (for some length scale L), for fre¬ 
quencies oj = n = jj. The mutual information of states 
and as well as the entropies of these states, 

approach different asymptotic values. 


shown in Fig. [3] In the limit of infinite acceleration the 
surviving correlations are given by I°°: 


TOO 

1 ^+ — 


8 In (2) 


(-In 


/531441\ 
V 256 ) 


(3 -f 2 ^ 2 ) In (3 -f 2 V 2 ) -h 

(3-2^2) 111(3-2^2)), 


TOO 

J.T,± — 


In 


FF In (2) VsVsy ■ 


(28a) 

(28b) 


As can be seen from Fig.[3l the entanglement entropies 
vanish for zero acceleration, as the mode is localized in 
region I. So there is no entanglement between modes 
in region I and modes in region II. As the acceleration 
increases, an acceleration horizon forms and the entan¬ 
glement entropy increases due to tracing out modes with 
support in the region behind the horizon. We see that the 
entanglement between modes in the accessible region and 
modes in the inaccessible region does not increase equally 
for all states, but depends on the particular ij^i. The 
mutual information of states 4’^^ and 4'^^ decreases 
with increasing acceleration and, in the infinite accel¬ 
eration limit, approaches distinct values (J^ ~ 0.4, 

I°°± ~ 0.6). Since I°° < I for all states, we conclude 

FF 

that also classical correlations become degraded with in¬ 
creasing acceleration. 

In this section we studied the degradation of quan¬ 
tum and classical correlations in fermion states that is 
caused by uniform acceleration. In the following section 
we study the entanglement in bosonic Bell states. 
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IV. ENTANGLEMENT OF UNIFORMLY 
ACCELERATED BOSON-BOSON STATES 


We continue by investigating the entanglement of Bell 
states of Unruh inodes uj and O of a massless uncharged 
scalar field, 


l^ls) =-^ (|0<^)c/|ln)c/ ± |la;)c/|0n)c/), (29a) 

l^ls) =-^ (|0<^)c/|0n)c/ ± |la;);7|ln)c/), (29b) 

I 


where uj, O are the frequencies and 0, 1 the occupation 
numbers of the Unruh modes. We consider the two modes 
Lj and n undergoing constant accelerations and ao, 
respectively. The acceleration parameters of the modes 
are denoted by r‘^ and r^. We write states (l29l) in the 
Rindler basis to obtain the infinite-dimensional density 

matrices ^ Then, to describe the system as it 

is seen by an observer confined to region /, we have to 
trace out modes that have their support in region II. 

As in Sec. cni to obtain the negativities W of states 
(1291) . we determine the partially transposed reduced den¬ 
sity matrices that are block diagonal and calculate 
the negative eigenvalues. More details can be found in 
Appendix[Bj The negativities of Bell states (1^ are given 
by the expressions 


N,,± =- 


1 1 


1 




Hz%r {z%y 

1 1 1 


^Z%Z% + \{n\ 




HZ%r (Zg) 


2 


LB, 


) + E^< 


(n.O) 



n—1 


m—1 


oo 

n—1 


(30) 

(31) 


where n% = (e^ — 1) ^ [ng = (e^u — 1) is the 
Bose-Einstein distribution with the Unruh temperatures 
Tu,/q = Z^^ is the bosonic partition function (fHl) . 

and 7 ,j,± is given by 

7$±g = 1 — (32) 

Ny}. , and give small corrections com- 

pared to the leading term and can be found in Appendix 
[bI The degradation of entanglement shows fundamen¬ 
tally different characteristics for the two Bell states 
and see Fig. [H While for entanglement van¬ 

ishes asymptotically, loses all its entanglement for 
finite acceleration. 

In case of for r‘^ = r^, only one of the blocks 

on the diagonal of the partially transposed reduced den¬ 
sity matrix admits negative eigenvalues. There is no en¬ 
tanglement generated in any sector, and only the sector 
where entanglement is initialized contributes to the neg¬ 
ativity . From now on, we refer to a sector as all 

B B 

elements of the reduced density matrix that contribute to 
one block of the block diagonal partially transposed re¬ 
duced density matrix. For r“ ^ there is entanglement 
created in all sectors. 

Note that the negativity vanishes asymptotically. The 
reason why 'kgg does not become nondistillable (and 
therefore does not become separable) for any finite accel¬ 
eration is that the occupation of state |00) is always zero; 
i.e., the Unruh effect does not drive the occupation of this 
state. This is due to the fact that the constituents of 


both contain one excitation and, thus, state |00) is not 
accessible. In consequence the matrix element |00)(00| 
in p^± is always zero and entanglement vanishes only 
asymptotically for infinite acceleration, as in this regime 
all occupation is shifted towards highly excited states. 
As we show in Appendix 111 the negativity of 

is of the form 7V^± = ^ 2 ^^ 771=0 A^?±™^ and it can be 

seen that each of the is bounded from above by 

The (partial) negativity can be 

^BB ^B B _ ^BB 

read off from m 


■^BB 


1 1 1 

^iZ%f {Z2f 





(33) 


where 7 *± is some kind of “cut-off function”. As all 

^BB 

]y{n m) bounded from above by , it follows that 

N^± vanishes for the same parameters as does. 

These parameters are characterized by = 1; i.e., as 

soon as this fraction of the population is excited to the 
first state above the vacuum, state loses its entan¬ 
glement. 

This can be understood in an intuitive picture consid¬ 
ering an effective state represented by an effective den- 

sity matrix p^yy{k) of the fcth sector. As we saw, the 
block diagonal nature of the partially transposed reduced 
density matrix leads to a negativity of the form 

R R 




E oo Ar(ra.O) Y^oo AI'I.U." 

n=l + Z^m=l 


^BB 


(0,m) 


So we 










FIG. 4: (Color online) Negativities for the maximally 
entangled boson states and where both ob¬ 

servers are accelerated, plotted versus the acceleration 
a = a^i = an, measured in units of j- (for some length 
scale L), for frequencies uj = fi = ^. The negativity of 
state (blue continuous line) vanishes for finite accel¬ 
erations, while the negativity of (red dashed line) 
vanishes asymptotically. 

introduce such that the negative eigenvalue of 

(^) provides Although there is not a 

strict symmetry between n and m in (ED, the effective 
description captures the essential features of the behav¬ 
ior of entanglement, as the vanishing of implies 

that is vanishing as well. Imagine that only one 

observer is accelerated ^ 0), while the other one is 
inertial (an = 0); then we can write an effective state 
/^e//(^) as 

Peffi^) =Ctfe(atcj)|Ofc)(l(/c -|- l)|-f 

+/?fe(aa;)|0(fc -|- 1))(0(A: -|-1)|-|- 
+^fe(a(^)|lfc)(lfc|-|- 

+7fe(a<.,;)|l(fc + l))(0fc|, (34) 

where 5k{a^) = 0 and we denoted \nQ) i by |nm). 

The coherences that are present in the initial state and 
are responsible for entanglement are quantified by Q:o(atj) 
and jo(O'oj)- The coefficients /3k{aa,) and Sk{ai^) have the 
physical interpretation of quantifying the occupation of 
the states |0(A:-|-1)) and \lk), respectively. Note that the 
Unruh effect drives the occupation of these states [in the 
present case only the occupation of |0(A: -I- 1))]. Initially, 
Pkiciuj) = Sk{aui) = 0. Now it is easy to see that, for 
fixed fc, (IMl) is always entangled for finite acceleration, 
but loses its entanglement for a^j —>- 00 , as in this limit 
ak^Pkilk -t 0 and thus Uk = Pk = h = Ik = 0. This 
explains why, for one accelerated observer (like in i), 
entanglement vanishes in the infinite acceleration limit 


but not for hnite accelerations. 

Moving to the general case of ^ 0, ao ^ 0, the 
effective density matrix of the fc-excitation sector is of 
the form 

ptff{k) =ak{a^,an)\0k){l{k + 1 )|-|- 

+PkiO'uj, an)|0(fc -l- 1))(0(A: -I- l)|-l- 
+5k{au,,aQ)\lk){lk\ + 

+ 7 fe(a<.j,aa)|l(fc-|-1))(0A:| (35) 

and the negativity vanishes for finite acceleration; i.e.. 
Oik = Pk = Sk = 7 fc 7 ^ 0 for a^, an < 00 . The equality 
between the strength of the coherences and the occupa¬ 
tion of states |0(fc -I-1)) and |lfc), i.e. ak = Pk = Sk = 7fe, 
is achieved due to the special structure of the reduced 
density matrix [cf. (IB2D ]. 

Let us have a look at the fc = 0 sector, where entan¬ 
glement is initialized. For vanishing acceleration there 
are the coherences | 00 )( 11 | and | 11 )( 00 | that are non¬ 
vanishing, while the states | 10 ) and | 01 ) are not occu¬ 
pied, i.e., Po = Sq = 0; this state is maximally entan¬ 
gled. By increasing the acceleration the symmetry of the 
term a^a^|nm)(nm| in the reduced density matrix p^± 
[cf. (IB2I) ] leads to an equal occupation of |10)(10| and 
|01)(01|. At the same time the coherences |00)(11| and 
|11)(00| are decreasing symmetrically. So for some hnite 
acceleration oq = ,do = <^o = 7 o 7 ^ 0 holds and entangle¬ 
ment vanishes. 

However, there is also entanglement creation in sectors 
of higher excitations fc > 0 that are initially unoccupied 
in the sense of ak = Pk = Sk = "fk = 0 . Although en¬ 
tanglement is initially increasing in these sectors due to 
acceleration, it is vanishing for the same acceleration as 
in the k = Q sector. This is due to the fact that, besides 
ap^app\nm){nm\, also a^a^|(n-|-l)(TO-|-l))((n-|-l)(m-|-l)| 
[cf. (IB2I1 ] drives the occupation of |lfc)(lfc| and there¬ 
fore compensates part of the loss of occupation of that 
state that is caused by the acceleration. This might be 
called “diagonal mixing” and is essential for achieving 
Ciklk = PkSk for hnite acceleration. This condition is sat- 
ished when rign^ = 1 holds. The fact that both modes 
“smear out” due to the acceleration enables Sk 7 ^ 0. This 
is the crucial point that enables the complete loss of en¬ 
tanglement for a hnite acceleration. 

Before moving on, we want to emphasize the depen¬ 
dence of the negativity on the energy of the inodes. The 
condition for entanglement is given by 

e~^ +e~"^ < 1 (36) 

and we see that entanglement is more persistent for 
higher frequency modes. We introduced the Unruh tem¬ 
perature T^/n = 2 ^^ ■ ff interesting to note that we 

can write condition (IMll . for the same acceleration for 
both modes, i.e., = uq = a, equivalently as 

w>riog (Z2)=-Fn, 


(37) 
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where T = ^ and Fq = — Tlog (Z^) is the Helmholtz 
free energy. Note that the same condition with uj and 
H interchanged also holds. So at least formally the 
Helmholtz free energy of one mode bounds the energy 
(frequency) of the other one. 

To summarize, by increasing the acceleration, i.e., by 
scanning through the families of states, entanglement de¬ 
creases for all bosonic Bell states. However, there is also 
entanglement created in sectors that have not been en¬ 
tangled initially. State loses all its entanglement for 
a finite value of the acceleration, whereas is entan¬ 
gled for all finite accelerations. This is due to the appear¬ 
ance of the function 7 *± that indicates the presence of 

^ BB 

a threshold, where the state becomes nonentangled. The 
reason for this behavior is the different occupation pat¬ 
terns of the constituents (structures) of the states we con¬ 
sidered here. The negativities for the states are plotted 
in Fig. [4] In contrast to the fermion case (cf. Sec. IHII) . 
all states lose their entanglement in the infinite accelera¬ 
tion limit. This is due to the infinite tower of excitations 
for bosonic modes that leads to a partition function Zb 
that grows unbounded. Intuitively speaking, the acceler¬ 
ation leads to a temperature that shifts the occupation to 
higher energy states and therefore the occupation of the 
lowest lying states approaches zero. Our findings provide 
evidence that the structure of the states plays an impor¬ 
tant role, as this decides about the set of states that are 
accessible. The “noise” introduced by the Unruh effect 
is state dependent. 

Next, after having addressed the bosonic case, we 
investigate the degradation of entanglement between a 
bosonic and a fermionic mode due to acceleration. 


V. ENTANGLEMENT OF UNIFORMLY 
ACCELERATED BOSON-FERMION STATES 

Using the same techniques as in Secs. [HT] and IIVI we 
study the degradation of entanglement in boson-fermion 
states. We start by considering the non-Bell states 

|Xi) (|0.)^|1^)+ + |Uc/|l^)0) , (38a) 

IA 2 ) (lU^|1^))) + IU()|1^)^) , (38b) 

where F labels the fermionic mode, w, O are the fre¬ 
quencies, and 0, 1 the occupation numbers of the Unruh 
modes. -I- and — refer to particles and antiparticles, re¬ 
spectively. The mode of frequency lj is bosonic while the 
mode of frequency Q is fermionic. The respective acceler¬ 
ation parameters are given by r = arctanh(e~~) for the 

TTfi 

bosonic and rj = arctan(e ) for the fermionic mode. 

Again we use the negativity (EH) as a measure of en¬ 
tanglement and obtain 


(39) 

(40) 


where Nf is the (universal) negativity that was found 
for maximally entangled fermions Nf = ^cos^{rf) = 
^{Zp)~^ [IH and Nb ^2 are given by 


Nb,i 

Nb,2 


11 

^ \Zb) 
_1 1 
"2^- 


(41) 

(42) 


These are the negativities in the case that only the bosons 
are accelerated f Appendix iBl). Details of the calculations, 
as well as the expression for Nn, can be found in Ap¬ 
pendix [cl 

Thus, the degradation of entanglement in states Xi 
and X 2 is quite similar to the behavior reported in 
Intuitively, what happens is the following. When acceler¬ 
ated the fermions get “rotated” and the bosons “smeared 
out”. Therefore, the fermions that are less affected by 
acceleration “mimic” the nonaccelerated bosons. On the 
level of the partially transposed reduced density matrices, 
we observe that the fermionic and the bosonic part fac¬ 
torize and thus the resulting negativity can be expressed 
in terms of negativities obtained from the cases of one 
accelerated observer. 

However, as we will see, this is not a generic feature 
and it is absent in cases of the boson-fermion Bell states 
and that are given by 

I^If) ± |0^)£/|ln)y) , (43a) 

I^If) ='^ (|0i.-j)c/|0a)(7 ± |l<.c;)£/|ln)F) ■ (43b) 


The negativities of states (031) are of the form N = 
Nn and again each of the Nn is bounded from above 
by Nq. Remember that Nq measures the negativity in 
the sector, where the entanglement is initialized. In the 
following we denote Nq by ,*± • The further 

for n ^ 0 can be obtained analytically (Appendix ICll . 
However, already with the expression for A^2± ,.± in 
hand, we are able to characterize N^± / 3 ,± . For states 

^bf/^bf 

'^If and we obtain 


where 






AtW = 


1 1 




=- 


1 1 


n—1 

(44) 

00 

n—1 

(45) 


(46) 

2 7$±^(nB,^F)- 

(47) 


Nx, =2NfNb,i, 
Nx, =2NfNb,2, 
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FIG. 5: (Color online) Negativities for the maximally 
entangled boson-fermion states X 2 , 
where both observers are accelerated, plotted against the 
acceleration a = = an, measured in units of j- (for 

some length scale L), for frequencies w = 17 = The 
degradation of entanglement occurs at different rates. 
The negativity of state (blue continuous line) van¬ 
ishes for finite accelerations, while the negativities of Xi 
(red dashed line), X 2 (purple, dotted), and (green 

continuous line) vanish asymptotically. 

Further, Zg, Zp are the partition functions (fldll . (fTSl) . 
ng = {e^ — 1)~^ is the Bose-Einstein distribution, 

n 

rip = (e^ -|-1)“^ is the Fermi-Dirac distribution, and 
the T^/q are the Unruh temperatures introduced by the 
acceleration. The function 7 *± is given by 

^ R IT 


(n) 

Details of the calculations and the expressions for 
and can be found in Appendix [Cl The ± 

fn) 

bound all the N / , , from above and therefore cap- 

ture the essential behavior of entanglement degradation. 
So in the following we restrict our discussion to these 
quantities and refer to them as negativity. 

In case of state the negativity is given by a prod¬ 
uct of the inverse partition functions for fermions and 
bosons. So, the negativity vanishes in the infinite accel¬ 
eration limit due to the unboundedness of the bosonic 
partition function. As for in Sec. lIVi the reason 

why N^l is positive definite for finite accelerations is 

'^SF 

that there are no contributions to the density matrix of 
the form |00)(00|, since these cannot be created by the 
Unruh effect for states 4'^. 

Moving to state 5*^^, we realize that, similarly to (1551) . 
the negativity of state vanishes for finite accelera¬ 
tions. Furthermore, as in (1331) . the threshold depends 


on the occupation numbers of the excited modes and 
n^. When the product ripTip equals 1, the negativity of 
state ^gp vanishes. So the threshold condition is of the 
same form as for ^gg, where it is given by ripTi^ = 1. 
However, why does the negativity of ^^p vanish while 
states Xi, X 2 , and '^^p are entangled for all finite accel¬ 
erations? First, we note that if either r or rj is vanishing, 
^gp is entangled for all finite accelerations. To find the 
answer for the generic case, we use an effective state de¬ 
scription, as in Sec. IIVI by 



where 

=ak{au,,an)\kO^){{k + l)l-^+|-|- 
+l3kio-uj, an)|(fc + 1)0’^)((fc + 1)0“^|-l- 
-\-Sk{au,, an)|^l"^~*')(fcl’^~*'| + 
+7fc(a<.j;ao)|(fc-I-1)1'^+)(A:0''^| (50) 

and 

=o^k{a^,an)\kl^~){ik + l)l^+l^"|-t 
+/3fc(fli.c;,on)|(fc + 1)1"^ )((^ + 1)1^ 1+ 
+Sk{a^,aQ)\kl^+l^~){kl^+l^~\ + 
+lk{auj,an)\{k + l)l^+l^“)(fcl^“ |, (51) 

where we denoted |ni^)/® |0n)j (g) |0n)7 by |n0^), |ntj)/(g) 
|ln)/(g)|0a)7 by |nl^+), |n^)/(g)|0a)|(g)|ln)7 by |nl^“), 
and |n„)/ (g) |la)| (g) |ln)7 by |nl^+l^“). Thus, if the 
fermionic mode is not accelerated, Sk = 0 and ^^p is 
entangled for all hnite accelerations. We see that the 
effective density matrix (1491) splits into two contributions. 
Further, it is easy to see that (l50l) and (jSTI) carry the same 
entanglement. Therefore, in the following, we restrict 

ourselves to Pg^^^(fc). 

Similarly to the bosonic state, ao = Po = lo = 
i5o ^ 0 for auj,an < 00 is achieved due to the spe¬ 
cial structure of the reduced density matrix p^± [cf. 
dClOp ]. For the further sectors, fc > 0, ak"fk = 
PkSk 7 ^ 0 is enabled. For vanishing acceleration there 
are the initially nonvanishing coherences |00^)(11^“'"| 
and |11'^+)(00^| that decrease with increasing acceler¬ 
ation. By increasing the acceleration some coherences 
are created {\k0^){{k + 1)1^''"| and |(A: -I- l)l'^+)(/c0'^|) 
and, further, the term i cos'^(r^)a^(|n0^)(n0^| -I- 
tan^(r/)|nl^+)(nl^+|) in (ICIOI) leads to an increasing 
occupation of |fcl^+)(A:l'^+| and \{k + l)0^)((fc -I- 1)0^|. 
In contrast to the bosonic state, this does not hap¬ 
pen symmetrically. However, at some point, when r and 
Tf fulfill ngTip = 1, there is an occupation of these two 
states such that akjk = PkSk 7 ^ 0 and entanglement van¬ 
ishes. 

Furthermore, as for the bosonic <I>^-state, the 
term icos^a^|(n -I- l)l^+)((n -|- 1)1^+|, as well as 
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i sin^(r/) cos^(r/)a^|nl^+)(nl^+| contribute to the oc¬ 
cupation of \{k -I- l) 0 ^)((fc + 1 ) 0^1 and therefore com¬ 
pensate part of the loss of occupation of that state that 
is caused by acceleration. Above we called this “diago¬ 
nal mixing”. This mixing enables akjk = Pk^k ^ 0 and 
therefore entanglement vanishes for finite accelerations. 

Again, we want to emphasize the dependence on the 
energy of the entangled modes. The condition for entan¬ 
glement can be written as 


< 1 - ( 52 ) 

This is equivalent to 7 $!^ > 0. As in Sec. HID we see 
that entanglement is more persistent for higher frequency 
modes. Furthermore, we note that, for equal accelera¬ 
tions = an), condition (15^ can be written in terms 
of the Helmholtz free energies as 


uj + n>T (log iZ%) - log (Z^)) = -F^ + Fo, (53) 

where denote the Helmholtz free energies. Com¬ 

paring that condition to (1371) , we see a huge similarity, as 
(l37l) can be written as w -I- H > —F^^ — Fq. If we would 
take these equations for more than just a nice rewriting, 
we could conjecture that the origin of the nonvanishing 
entanglement for fermions is given by the fact that the 
condition for entanglement is lo + n > F^^ + Fq, where 
Fui + Fq < 0 , and thus it is trivially fulfilled for all accel¬ 
erations. 

To summarize, the negativities of states Xi and X 2 
factorize and we observe a product structure similar to 
the one obtained in Sec. cni where the total negativity 
is the product of the fermion and the boson contribu¬ 
tions. That is due to the structure of the fermion mode 
n. These families of states are entangled for all finite 
accelerations. The negativities are given by the product 
of inverse bosonic and fermionic partition functions and 
therefore vanish in the limit of infinite acceleration. In 
case of the and states (‘f’^j’ ^bf) 'i®®® 
not hold any more and we observe a behavior that is sim¬ 
ilar to the one we obtained for the and states in 
Sec. IIVI Again, state $ 3 ^ loses all its entanglement for 
finite accelerations, while state ' 1 '^^, is entangled for all 
finite accelerations; see Fig. [5l The different behavior is 
due to the different structures of the states, since only in 
the case of state diagonal mixing is enabled. 

So, we have seen that, for fermion-fermion, boson- 
boson, and boson-fermion Bell states, the degradation of 
entanglement does not depend on, for example, whether 
the state is a singlet or a triplet, but on the structure of 
the particular state. It is the structure of the state that 
determines the fading of its entanglement. 

In the following section, we summarize our findings and 
discuss the role of particle statistics in the degradation 
of entanglement. 


VI. ENTANGLEMENT DEGRADATION AND 
THE ROLE OF PARTICLE STATISTICS 


In this section we discuss the mechanisms behind en¬ 
tanglement degradation and the role of particle statis¬ 
tics therein. Above we discussed the fermion-fermion 
Bell states (fTCl) . the boson-boson Bell states (1^ and the 
boson-fermion Bell states (H51) . Using the expressions for 
the negativities (1^ . (fTOl) . (1501) . and (ED, as well as (HU) 
and dUl), we can write the negativities of all Bell states 
in a compact form 



I I I 


(54) 


where S^y = denotes the entangled state. A, 

Y encode the statistics of the fields (fermionic, bosonic), 
and x,y are equal to 1 for fermions {X,Y = F) and 
equal to 2 for bosons {X,Y = B). The functions 7 g± 
are given by 


=\J - nB, (55a) 

=1, (55b) 

=^Z%Z^ + {Z^Z^fnl - Z^Z^fiF (55c) 


for the '^xY states, where Ub/f = ^{^BjF +^® 
the average occupation number, and 




=l-n‘ 


B'‘‘Bi 


^ BF 



(56a) 

(56b) 




= 1 , 


(56c) 


for the ^xY states. Equation (l54)l gives the negativities 
in the fc = 0 sector, where entanglement is initiated. For 
some states there is entanglement dynamically created in 
other sectors but these are always bounded from above 
by (l54l) : see Fig. [H Therefore, these negativities capture 
the main features of entanglement degradation and we 
restrict our attention to these. 

Physically speaking, after fixing the frequencies w and 
n, as seen by the accelerated observers, Eq. (l54ll gives the 
negativity of the two-parameter family of states Sxy 
That is, for each choice of the pair (atj,an), Eq. (15^ 
gives the negativity of the particular state S^y f^at is 
characterized by (w,fl,a^,an), when this state is seen 
by accelerated observers of accelerations a^j and an, re¬ 
spectively. For this reason, and also because the Unruh 
modes that are considered are global modes, the setting 
should be considered as a toy model that captures the 
essential features of entanglement degradation. 

There are essentially two factors determining the fad¬ 
ing of entanglement. The first one is given by the set of 
states (as above, by states we mean diagonal elements of 
the density matrix) that become available when a state is 
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(a) State for accelera- (b) State ^gg for accelera¬ 


tion flij = an = 50-i 


tion au, = an — 100-^ 




frequency co 


frequency co 


(c) State $ 3 ^ for accelera¬ 
tion Ou; = an = 50i 


(d) State accelera¬ 

tion aui = an = 100 i 


FIG. 6: (Color online) Energy dependence of the entan¬ 
glement in Bell states ^bb/bf- States of modes of fre¬ 
quencies oj and n (in units of for some length scale 
L) contained in the red region show no entanglement, 
while states of higher energies remain entangled (white 
region). The plots show states for acceleration 

aui = an = 50;^ ((a), (c)) and a^, = an = 100-^ ((b), 
(d)). State is entangled for all frequencies and ac¬ 
celerations and therefore it is not shown. The zeros of the 
functions 7 ,j,± and 7 *± define the border between the 
regions. The asymmetry in (c) and (d) is due to the fact 
that fermions are “more resistant” towards the effects of 
acceleration. 


accelerated. This set depends heavily on the structure of 
the state, as, for example, for the state |00)(00| never 
becomes accessible, but also on the statistics that does 
not allow for two- or more-particle states for fermions. 
The second determining factor is whether the population 
of states |00)(00|, |01)(01|, |10)(10|, and |11)(11| can be 
transferred completely to higher excited states like, for 
example, |21)(21|. If that is possible also the coherences 
|00)(11| and |01)(10| vanish with increasing acceleration. 
Both factors depend heavily on the statistics of the un¬ 
derlying field. 


For illustrating reasons, consider the density matrix p. 


o 

o 

0 

0 

Pooil \ 

0 

Poioi 

Poiio 

0 

0 

PlOOl 

PlOlO 

0 

V PllOO 

0 

0 

$ 

Pnii / 


written in the basis {|00), |01), |10), |11)}, where all 
{i,j,k,l G {0,1}) are zero for states and vice versa. 
Tracing out anti-particles, (EZD is the full density ma¬ 
trix for fermions and acceleration decreases the diagonal 
elements (Poooo^ PfiiF Poioi’ Pioio)’ as well as the co¬ 
herences (pfoii; Pnoo; Pwoi’ Poiio) to a finite value for 
both <I>^ and states. Furthermore, the “squared” co¬ 
herences PooiiPfioo and PipoiPmio always dominate the 
product of the occupations of states Poooo’ pfiii and 
PmoiJ Pwioy respectively. That is why entanglement de¬ 
creases but does not vanish. 

For bosons, in contrast, (l57ll is only the initially non¬ 
vanishing part of the infinite-dimensional density matrix. 
The reason why states lose their entanglement (cf. 
Fig. [S]), while are entangled for all finite accelera¬ 
tions, is given by the set of available states (cf. Secs. ITVl 
and 0. The reason for the asymptotic vanishing of the 
negativity for all boson-boson and boson-fermion states 
is given by the second determining factor (see above). All 
initially nonvanishing elements of p, i.e., the and the 
pfjkF approach zero in the infinite acceleration limit, as 
matrix elements of higher excitation states are increas¬ 
ing. That is due to bosonic statistics. Naively, fermions 
get “rotated” and bosons “smeared” towards higher ex¬ 
cited sectors due to acceleration. Furthermore, note the 
asymmetry in Fig. [S] that is due to the fact that fermions 
are “more resistant” towards the effects of acceleration. 

This explanation based on the two determining factors 
captures the role of particle statistics. Particle statistics 
also is reflected in the negativities (I5H) that are written 
in terms of the partition functions and the occupation 
numbers and therefore make the “effects of statistics” 
apparent. 

In this work we studied maximally entangled states 
like, for example, 

|^'„) =sin(a)|0<^)(7|ln)(7-I-cos(a)|l^)j/|0a)j/, (58a) 

|$„) =sin(a)|0<^)(7|0n)(7-I-cos(a)|l^)j/|la)j/, (58b) 

where we chose a to be 7r/4. Differing choices of a lead 
to less entangled states, as the negativity = 

sin(Q!) cos(a) is maximized for a = 7r/4. Given the mech¬ 
anisms that lead to entanglement degradation that we 
outlined above, it is evident that states (l58l) behave qual¬ 
itatively the same for generic a as they do for a = 7r/4, 
i.e., in the maximally entangled case. In the case of more 
general mixed states, we expect that, depending on the 
state, one observes that the degradation of entanglement 
shows characteristics that are best described by a mixture 
of the characteristics of the degradations in the case of 
and However, it seems reasonable to expect that 









13 


a random mixed state will lose its entanglement for finite 
acceleration with high probability, as in the fermionic 
case studied in (4ll |. 

In some sense we can think of the functions 7 [ (1551) and 
(l56l) ] as deformations of a (universal) particle-statistics- 
dependent negativity Mu , 


Mu = 


1 1 


1 


2 {z-^r (z^) 


n\y ’ 


(59) 


that only depends on the partition functions that are 
characteristic for the particle statistics of the field. Then 
the particular structure of the Bell state, as well as the 
particle statistics, set the 7 that we might call structure 
functions and denote them by jstructure- As we saw, 
these depend heavily on the set of states whose occupa¬ 
tion is driven by the Unruh effect. So finally we can write 
the entanglement of a Bell state (negativity Nstate) in the 
sector where entanglement is initialized as 


The entanglement of states $^5 'i’BF com¬ 
pletely vanishes at a finite distance from the horizon that 
is large compared to the Planck length Lp {d ~ 0.01i?s, 
where Rs is the Schwarzschild radius), whereas states 
are entangled for any finite distance from the black 
hole. 

Thus, we observed that entanglement, an important 
resource for quantum information tasks, gets degraded 
very differently for differing Bell states, i.e., the degra¬ 
dation of entanglement is state dependent. Our findings 
imply that there are particular states that remain entan¬ 
gled as seen by an observer that is uniformly accelerated 
or equivalently is stationary close to a black hole, while, 
for other choices of the state, there is no entanglement re¬ 
maining. This implies that the gravitational degradation 
of entanglement depends on the structure of the state. 


VIII. CONCLUSIONS 


Mstate — ^structure Mu • (^ 0 ) 

We want to close this section by giving some comments 
on Eq. O- First, choosing the partitioning of Ngtate in 
Istructure and Mu is not unique and there are possible 
partitions different from (1541) . Nevertheless, writing the 
negativity in form (1601) makes the importance of the par¬ 
ticular structure of the state manifest. Moreover, Eq. 
( 1601 ) makes it possible to clearly identify the two deter¬ 
mining factors: The first one sets the function structure, 
while the second one determines Mu- Further, we note 
that a slightly varied form of (1501) also holds for states Xi 
and X 2 - We expect that there are slight modifications 
when there are different particles involved, like particles 
carrying spin. Finally, it would be interesting to figure 
out whether the negativities of non-maximally entangled 
mixed states could also be captured in an expression sim¬ 
ilar to (l60|) . 

In the following section, we briefly point out possible 
implications of the above findings for particles in Bell 
states close to the black hole horizon. 


VII. DEGRADATION OF ENTANGLEMENT IN 
THE VICINITY OF A BLACK HOLE 

The framework we used in this work also applies to 
the spacetime close to a black hole, as outlined in Ap¬ 
pendix [D] (see also [l^). However, given the caveats in 
the interpretation of the states, we described in Sec. m 
the following discussion aims at giving a qualitative idea 
about entanglement degradation near black holes. 

One point that can be inferred is that entanglement 
gets degraded in the vicinity of the black hole horizon. 
Further, all states that involve bosons lose their entangle¬ 
ment in the limit of reaching the horizon. This is in con¬ 
trast to the fermion-fermion states, where entanglement 
never vanishes. Further, there are crucial differences be¬ 
tween the degradation of entanglement for states and 


In this work, we studied families of two uniformly ac¬ 
celerated maximally entangled Unruh modes in the gen¬ 
eral case of different accelerations and analyzed their en¬ 
tanglement, measured by the negativity. Therefore, we 
considered states containing two fermionic modes, two 
bosonic modes, as well as states of one bosonic and one 
fermionic mode. Special emphasis was given to the com¬ 
parison of Bell states and 4'=*=. Although the Unruh 
modes we used do not have a simple physical interpre¬ 
tation, our studies provide insight into the mechanisms 
that lead to the degradation of entanglement due to ac¬ 
celeration. 

We found that, in contrast to the other cases, purely 
fermionic families of Bell states are entangled for all ac¬ 
celerations. Still, the entanglement of state 4^^ degrades 
faster with acceleration than the entanglement of state 
Interestingly, it is only for state 41^ that both ac¬ 
celerated modes give rise to a contribution to the same 
diagonal element of the reduced density matrix that is 
relevant for entanglement. We suspect that this special 
feature of 4/^ is responsible for the different behavior 
regarding entanglement degradation. Furthermore, we 
found that also classical correlations are partially lost 
due to acceleration. 

In the purely bosonic case, as well as in the boson- 
fermion case, state 4'=*= remains entangled for all finite 
accelerations, and entanglement vanishes asymptotically 
in the limit of infinite accelerations. In contrast, state 
loses its entanglement for some finite acceleration. This 
is manifest in the presence of a “cut-off function” 7 $± in 
the expression for the negativity. So we found that the 
type of Bell state (i.e., being 4>^ or 4>=*=) crucially affects 
the robustness of its entanglement against acceleration. 
Furthermore, we obtained that the reason for the occur¬ 
rence of this phenomenon is originated in the particular 
occupation patterns of the constituents (the “structure”) 
of the state, which determine which excitations can be 
driven by the Unruh effect. 
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Applying an effective state picture, we were able to 
explain this crucial difference between both types of 
states. State is entangled for all finite accelerations 
as the Unruh effect does not drive the occupation of state 
|00)(00|, and this state is naturally absent in the density 
matrix of 4'^. Entanglement vanishes only asymptoti¬ 
cally for infinite acceleration, as in this regime all occu¬ 
pation is shifted towards highly excited states. For state 
4)^ things are different. Essentially what happens is the 
following. For vanishing acceleration, there exist the co¬ 
herences |00)(11| and |11)(00| that are responsible for 
the entanglement, while the occupation of |10)(10| and 
|01)(01| is vanishing. When acceleration is increasing, 
the coherences are decreasing, while at the same time 
the occupation of |10)(10| and |01)(01| is driven (sym¬ 
metrically) by the Unruh effect by creating one excita¬ 
tion in |00)(00| . Thus, for the value of the acceleration 
for which the condition for entanglement [cf. (1361) . (15^ ] 
is violated, entanglement vanishes. Hence, we traced the 
difference in the behavior regarding entanglement degra¬ 
dation back to the set of accessible states and the sym¬ 
metry in the distribution of probability among them. It 
seems that diagonal mixing, as we coined it above, is re¬ 
quired to achieve sufficient uniformity in the occupation 
of the states. 

Further, we found that there are two factors that de¬ 
termine the fading of entanglement. The first one, given 
by the set of states that become accessible due to the 
Unruh effect, is heavily influenced by the structure of 
the state. Said factor determines whether a state loses 
its entanglement for finite accelerations. The second fac¬ 
tor is more closely related to the particle statistics of the 
modes that constitute the Bell state. It is whether higher 
excitation states become accessible due to acceleration. 
That is the case for bosonic modes, and thus Bell states 
in which these modes are involved are nonentangled in 
the infinite acceleration limit, whereas purely fermionic 
Bell states are always entangled. Remarkably, we found 
that the negativities of the boson-boson, boson-fermion, 
and fermion-fermion Bell states can be expressed in the 
same form (l54|) . 



1 I I 

2 (z^)^ 


(61) 


where the are the partition functions (of a har¬ 

monic oscillator or two-level system with energy gap 
and the 7g± are functions determined by the first 
factor, we introduced above. 

Furthermore, we discussed possible effects of hovering 
over a black hole on entangled states of two Unruh modes. 

In summary, our studies reveal the mechanisms that 
cause the behavior of entanglement in accelerated frames 
to depend heavily on the particular occupation patterns 
of the constituents of the entangled state. 
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Appendix A: Eermion-fermion states 


1. Calculation of negativities 

The negativity for a composite system (we denote the 
subsystems by A and B) described by a density matrix 
P = Pab is given by the sum of the absolute values of the 
negative eigenvalues of the partially transposed density 
matrix p^Jg^ 


3 

where the A^ ’s are the eigenvalues of p^g and for a den¬ 
sity matrix pab = Y.ki 7 nnPkimn\k){l\®\m){n\ the partial 
transposed is given by p^Jg = Y.kimn Pki mn |fc)(/|(g)|n)(m|. 

The calculations of the negativities of the fermion 
states {'0i} = have to be carried out with 

care and the braided tensor product has to be taken into 
account [43| . As above, when we introduced the Unruh 
modes, we chose the ordering \ijkl)i:j = \ia,)^ 0 <8 

® The density matrices are obtained as 

Pj jj = To obtain the reduced density matri¬ 
ces p, = Trji {p%i') (i e trace out 

modes supported in region II and take care of the oper¬ 
ator ordering. Finally, we partially transpose the reduced 
density matrices and identify the blocks of p^'^ that ad¬ 
mit negative eigenvalues. 

In the case that the observers are able to detect parti- 
cles as well as anti-particles, the relevant blocks of 

P^I+ are given by 

^ FF 

(A2) 

62^^ =tan2(r;)cfEi, (A3) 

63"'-=tan2(r^)ctV, (A4) 

64 =tan^(r^) tan^(rj)cE^!^i; (^ 5 ) 
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where 


the negativity as 


=^cos^(r;)cos"(ry) 


icos{rJ) cos{r^) 


ecos{rJ) cos{r^) \ 
sin^(rj) +sin^(rj) J ' 

(A6) 


iV^i^ =(1 + tan^(rj) + tan^(r5^) + tan^(rj) tan^(rj ))x 

xN{c%n 

=i cos^(rj) cos^(r^). (A14) 


The negativity N^+ is sum of the absolute values of 
the negative eigenvalues of . Thus, using that 

and the fact that the can 

be obtained from the by the replacement e —^ —e, 

we can write 

iV,j,± =(1 + tan^(r|^) + tan^(rj ) + tan^(rj) tan^(rj ))x 

xN{c%n 

= \i- {sin^{r‘f)+sm'^{r^)) + 



(A7) 


As expected, assuming that only particles can be de¬ 
tected by the observers, i.e., tracing out anti-particles the 
negativity does not change. In this case, there is only one 
block contributing to the negativity that is given by 


Assuming that the observers can only detect particles, 
does not affect the negativity. Tracing over anti-particles 
in leads to 

1 / sin^(rj) cos^(r-^) cos(rj) cos(r^) \ 

2 \ cos(r^) cos(r^) sin^(rj ) cos^(r|^) ) 

with the negative eigenvalue —fV^± . 


2. Entanglement in different sectors 

We briefly discuss the distribution of entanglement 
in states '0^, when the acceleration is non-vanishing. 
While entanglement in the sector where it was initiated 
is decreasing with increasing acceleration, there is en¬ 
tanglement created in previously nonentangled sectors. 
Schematically, we can write fV,i,± as 


1 / 0 cos(r|^) cos(rj) \ 

2 \ cos(rj) cos(rj ) sin^(rj)-|-sin^(r-^) y ’ 


(A8) 


which has the negative eigenvalue —N^± and thus en- 
tanglement remains unchanged. 

Considering state in the case that the observers 
are able to detect particles as well as anti-particles, the 

relevant blocks of are calculated as above and 

^FF 

are given by 


bf-- =tan2(r^)ct%i, 

bf-- =tan2(r5?)ctT, 

64^^ =tan2(r^)tan2(r;;^)cf|^4. 


(A9) 

(AlO) 

(All) 

(A12) 


where 

=i cos^(rf) cos^(r9)x 

2 J j 

/ sin^(rj) cos^(rj) e cos(rj) cos(r^) \ 
y ecos(rj) cos(rj ) sin^(rj ) cos^(ry) y 

(A13) 


We notice that iV(cfA0') = and that again 

—>■ is induced by e —>■ —e. So we can write 


=iV(0,0|l+,l+) 

-f7V(l-,0|l+l-,l+) 

+iv(o,i-|i+,i+i-) 

+fV(l-,l-|l+l-,l+l-), (A16) 

where 

lV(0,0|l+,l+) = lV(ctT), (A17) 

A^(l-,0|1+1-,1+) =tan2(r}^)7V(cfif), (A18) 

iV(0,1“ 11+, 1+1“) = tan2(r5*)A^(cfJf), (A19) 

A^(l-,1- |1+1-,1+1-) =tan2(r}?)tan2(r^)7V(cf|f) 

(A20) 

are the negativities of the different sectors. 
A^(0,0|1''",1''") is the negativity of the sector, 
where the entanglement is initialized, i.e., entangle¬ 
ment between |0„)| ® |0^)7 (g) |0n)j ® |0n)7 and 
|l(.j)/ (g) |0t„)7 (g) |ln)j 0 |0n)7- There are three sectors 
in which entanglement is created due to acceleration 
(IAl8ll - (IA20ll . The negativities are plotted in Fig. [7] 
and we see that the entanglement is distributed equally 
between the sectors in the infinite acceleration limit. 

For state we see again that entanglement is cre¬ 
ated in some sectors and we can write N.^,± schemati- 
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cally as 

=7V(1+,0|0,1+) 
+Af(l+1-,0|1-,1+) 
+Af(l+,1- |0,1+1-) 


+fV(l+l“,l“|l“,l+l“), 

(A21) 

where 


iV(i+,0|0,l+) = Ar(ctV), 

(A22) 

iV(l+l“, 0 1 1“, 1+) = tan2(r7)iV(ctT), 

(A23) 

Ar(l+, 1“ 1 0,1+1“) = tan2(r^)iV(cfif), 

(A24) 

A(l+l“,l“ 11“,1+1“) = taii^(r‘^)tan^(rf 

(A25) 


are the negativities of the different sectors. 
fV(l+,0|0,l+) is the negativity of the sector, where 
the entanglement is initialized. When the accelera¬ 
tion increases from zero the negativity 0 | 0,1’*') 

decreases while (|A23ll - (IA25l) increase. In the infinite 
acceleration limit all sectors are equally entangled; see 
Fig. [3 Comparing states and we note that 

the “redistribution” of entanglement is the same for 
small accelerations, but differs more and more with 
increasing acceleration, until finally different limiting 
values are approached. 

We plot the entanglement in the sectors that show 
creation of entanglement in Fig. [T] This effect is due 
to a symmetric production of particle-antiparticle pairs. 
We can, for example, write state schematically as 
0, 0|1“'',1+. Then, for nonzero acceleration, the Unruh 
effect, that in some sense can be seen as pair produc¬ 
tion, populates states 0,1“ and 1+, 1+1“ and thus cre¬ 
ates entanglement [Af(0,1“ 11+, 1+1“) ^ 0] in the sec¬ 
tor 0,1“|1+,1+1“. In the same maimer there is en¬ 
tanglement created in the sectors 1“,0|1+1“,1+ and 
1“, 1“ I 1+1“, 1+1“; see Fig. [3 However, in general, not 
all states show this behavior. One state that does not 
show this feature is {\1^)u\^q)u + \l!^)u\^n)u)■ The 
reason why there is no entanglement generated is that for 
this state symmetric pair production would require a vi¬ 
olation of the Pauli principle and therefore is forbidden. 

Appendix B: Boson-boson states 

We calculate the negativitiy (lAll) for the bosonic Bell 
states 

=-^ (|0<.,;)(7|ln);7 ± |li..;)(7|0n)(7), (Bla) 

I^Ib) =-^ (|0<.,;)(7|0n);7 ± |li..;)(7|ln)(7), (Bib) 

where w, are the frequencies and 0, 1 the occupation 
numbers of the Unruh modes. The two modes cj and U 



FIG. 7: (Color online) Entanglement distribution for 
states (red dotted line) and (blue continuous 
line) plotted versus the dimensionless acceleration 
parameter r = = r^. The infinite acceleration 

limit corresponds to r = ^ . The middle pair of 
solid curves is twofold degenerate due to the sym¬ 
metry of the states. With increasing acceleration, 
entanglement is created in previously separable sec¬ 
tors. From top to bottom the curves correspond 
as follows: red, A^(0,0 11+, 1+), A(l“, 0 | 1+1“, 1+) 
and A^(0,l“|l+,1+1“), A(l“,l“ 11+1“,1+1“); 

blue, Ar(l+,0|0,1+), fV(l+l“,0|l“,l+) and 

iV(l+, 1“ I 0,1+1“), fV(l+l“, 1“ 11“, 1+1“). 


undergo constant accelerations a^i and ap. The accelera¬ 
tion parameters of the modes are denoted by and r^, 
respectively. 

We start with state and for concreteness we con¬ 
sider ^pp that shows the same entanglement as ^pp- 
So in the following we always denote ^pp by ^^p and 
similarly for the '!'+ states. The density matrix p^± 
that we obtain after tracing out region II is given by 

m,n 

+ l + l)(w + l))((n + l){m + 1)1 

m,n 

+ 2 ^ CLm0.n(im0.n\ nm){{n + l)(m -I- 1)| -I- h.c., 

m,n 

(B2) 

where |nm) = \nui)i 0 |ma)/ and a„ = 
an{r‘^) = tanh"(r‘^) cosh“^(r‘^), a„ = a„(r“) = 

tanh"(r‘^) cosh“^(r‘^)-\/rrTT. Accordingly, = 

an=m(r^) and am = an=m(r^)- The negativity can be 
obtained as a sum over the negativities for different 
values of n,m (JV^± = mThis can be 

seen from the block diagonal structure of the partially 
transposed density matrix. The part of the partially 
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transposed density matrix that contributes to 

I 


and the negativity are given by 

^ R R 


ltanh^"(r‘^)tanh^’"(r-^) 
2 cosh^ {r ‘^) cosh^ {r ^) 


1 + 


(m+l)n 

sinh^ (r‘^) sinh^ (r^^) 

y/{m+l){n+l) 
cosh(r‘^) cosh(r^^) 


^ tanh^(r^) 


1 + 


y/ (m+l)(ri+l) 
cosh(r‘^) cosh(r‘“) 
(n+l)m 


sinh^(r‘*^) sinh^(r“) 


) tanh^(r“) 


(B3) 


and 




-E^fe7 = E<r+E^. 


r(0,m) 


E 


BB 
2n I 


tanh^"(r‘^) 


=iV*'± , / , 

‘*’bb 4 cosh^ (r ^) cosh^) 


(tanh^(r“) + tanh^(r^) + 


sinh^(r‘^) cosh^(r^) 


(— cosh"^(r‘^) sinh"^(r*^) tanh"^(r”) + sinm(r‘^) sinlv(r*^)x 


sinh^(r“) sinh^(r^) cosh(r“) cosh(r'^) 4 

x(^ + - cosh^(r‘^) sinh^(r^) tanh^(r^) + 1) + ^ sinh^(r“) cosh^(r“) sinh^(r2) tanh^(r^)+ 

1 

+ sinh®(r“)(- tanh^(r‘^) sinh^(r^) cosh^(r^) — - sinh®(r^)))2) + (w O 11; n —>• m), (B4) 


where we used that ^ 0 only for either n = 0 or 

m = 0oTn = m = 0. It can be seen that each of the 

is bounded from above by N^2 = that de- 

‘J-BB «’Ib ‘^bb . . . 

scribes the entanglement between the modes we initially 
start with. For non-vanishing acceleration there is entan¬ 
glement created between higher modes, i.e., 7"’™^ ^ 0, 

but this will be a small contribution compared to ; 

^BB 

see Fig. [Sal 

To obtain , we set n and m to zero in (ESI. The 
negativity can be written in the form 

^BB 


iv7 =2iVo.^iVo,o7$±, 


(n%, n§) , (B5) 


where Uq = — 1)“^ is the Bose-Einstein distri¬ 

bution and Nq i^/q denotes the negativities if only mode 
w/fl is accelerated. These can be obtained from (IB4I1 by 
setting the acceleration parameter to zero, 


E± 

■^BB 




(B6) 


where is the bosonic partition function (1141) . Fur¬ 

ther, 7^± is given by 


7,(,±^ = l-n^ng. (B7) 

Now we move to state where the relevant part of 
the partially transposed reduced density matrix is given 
by the following expression: 


1 tanh^"(r“) tanh^’"(r^) 

2 cosh^ (r‘^) cosh^ (r^) 


n V (m-l-l)(n-|-l) 

sinh^ {r^) cosh(r^) cosh(r^7 

V O'^v^A^+l) ( I .. N tanh^(r") / . n tanh^Cr^^) 

cosh(r‘^) cosh(r^7 ^ ' cosh^(r^) ' cosh^(r‘^) 


sinh^(r^7 

y/{rn+l){n+l) 


(B8) 


Contrary to (ESIi only the eigenvalues of the block {n = 
0,m = 0) can be negative [Fig. |8b] for = r^. In 
this case, the sum of the partial negativities N^± = 

^ R R 


J2n m collapses to N^± = and we find 

’ '^BB '^BB 
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(a) Negativities 





FIG. 8: (Color online) Negativities, where both observers are accelerated, plotted against the acceleration a = a^j = ao, 
measured in units of ^ (for some length scale L), for frequencies oj = ^ While for states and 

there is entanglement created in initially nonentangled sectors (a), (c), and (d), there is no entanglement production 
(with this special choice of the acceleration parameter r = r‘^ = r^) for (b). Therefore in the generic case there 

is entanglement generated in initially nonentangled sectors. 


In the case ^ we assume without loss of generality 

I 


> r^. Then the blocks for m = 0 admit negative 
eigenvalues and we find the negativity 


N..± =N„p 


r(0) 


E 


tanh^"(r‘^) 


(—n coth^(r‘^) cosh^(r“) — (n + 1) sinh^(r“) — sinh^(r“)+ 


2/ n 


2^'„n^ 


2 /^u)\ 


^ 4 cosh"^ {r‘^ ) cosh^ (r^) 

H-2-(f^ cosh(2r“) cosh(2r^) + — + sinh^(r‘^) sinh^(r*^) + sinh^(r‘^)') + 

sinh (r"^) ^2 2 / 

— - sinh^(2r‘^) sinh^(2r^) (n{n + 1) sinh^(r*^) — sinh^(r‘^) sinh^(r^)))2)j 


(BIO) 


where is given by (IB9D . 


Appendix C: Boson-fermion states 

We calculate the negativities for maximally entangled 
states of a bosonic mode entangled with a fermionic one. 
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To start, we consider the states 


is, 


|Xi) (|0.)j/|l^)^ + \Uu\lE)u) > (Cla) 

1X2) {\Utj\l^)u + \Ud\^n)i) , (Clb) 

where F labels the fermionic mode, w, 12 are the fre¬ 
quencies and 0, 1 the occupation numbers of the Unruh 

modes. -I- and — refer to particles and antiparticles, re¬ 
spectively. The mode of frequency w is bosonic while the 
mode of frequency 12 is fermionic. The respective acceler¬ 
ation parameters are given by r = arctanh(e~^) for the 

TTfi 

bosonic and Vf = arctan(e ) for the fermionic mode. 
As in Sec. El we carefully take into account the operator 
ordering for fermions, especially when we are performing 
partial traces. The relevant part of the partially trans¬ 
posed reduced density matrices for these states can be 
computed to be 


1 . 1 2n / \ / / n+1 

1 2 / \ (^) I sinh^ (r) y cosh^ (r) 

2 "^°" cosher) tanh2(r) 


for state Xi and 

1 Aanh2'"+2-(r) 

- cos (r-f)-=— 

2 cosh®(r-) 


n 


V^(m-h l)(n-t-1) 
yq?7i + l)(n+l) m 


(C2) 


(C3) 


Nx, =2NfNb,u (C4) 

Nx, =2NfNb,2. (C5) 

where Nf is the negativity Nf = 2cos^(r^) = 
and Nb^i, Nb ,2 are given by 


11 

^ \'^b) n=l 


^ cx; 

Nb,2 =- - -3 ^ e 

\^B) n,m—Q 


27r(n + Tn) 


1 1 
2^’ 


(C6) 

(C7) 


i.e., the ones we obtain when we only accelerated the 
bosons. The Nn in (IC6II can be obtained by setting = 
0 in (IB3I) and are given by 


iv. 


2cosh^('r) 2sinh^(r'l 2 


-I- - tanh^(r)-|- 


^ + \ sinh^(r) tanh^(r) -|- 1 


4sinh^(r) 


cosh^(r) 


)■ 

(C8) 

Next we move to Bell states and that are 
given by 


I^bf) ='^ (|1c.j)(7|0q)c/ ± |0^)(7|1q)^) , (C9a) 

I^bf) ='^ (|0<.j)(7|0q)b ± I!(..;) (71 In )b) ■ 


for state X2. We observe that the fermionic and the 
bosonic part factorize and so we find that the resulting 
negativity can be expressed in terms of negativities ob¬ 
tained from the cases of one accelerated observer. That 


The negativity of state is calculated in the following. 
For concreteness we carry the calculations out for state 
®BF- After obtaining the reduced density matrix p^+ 
by tracing out region //, 


P<i.+p = 2 ^^"' iYl )QnI»)(»I ® (|00)(00| -f tan2(r/)|10)(10|) -f ^d^|(n-f l))((B-f 1)1 ® |10)(10| 

n n 

+ ^ cos(r/)a„d„|B)((n -f 1)| (g) |00)(10| -h h.c ■nondiag. } 
n 

irf) { cos^{rf)al\n){n\ g (|01)(01| -f tan2(r/)|ll)(ll|) + ^ |(n -f 1 ))((b -f 1)| g |11)(11|-H 

n n 

'y ^ COs(rj )n^O,yj In) ((n-t- 1)| g |01)(11| -t- h.c. nondiag. } j (^70) 

n 


where a„ = a„(r) = tanh"(r) cosh~^(r), an = d„(r) = 
tanh"(r) cosh~^(r)-\/n -|- 1, and the notation \ij) = 
Ka)/" g bn)/■ Then after partial transposition, the rel¬ 
evant part of the reduced partially transposed density 


matrix is of the form 


V 0 tan2(r/)cfjfi / 


(Cll) 
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where 


$+ 1 2/ .tanh^”(r) 


cos2(r-/)tanh2(r) e cos(r/)^/^^^ 
ecos(r/)y^^g^ ^^+sin\rf) 


■ (C12) 


Already at this stage we see that the fermionic contri¬ 
bution does not simply “factor out” as it was the case 
for states Xi and A2. Due to the block diagonal form 
of (ICllI) . the negativity A^$+p that equals can be 

written in the form 


; see Fig. [Sc] To obtain we have to calculate 

^BF ^BF 

^i± 

^BF 


=2NfNo,^e^ 

^BF 


n^n% - npn% ) , (C14) 


where n% = — 1) ^ is the Bose-Einstein distribution, 

Nq^cj is given by (IB6I) . and rip = (eAT -(-l)“i is the Fermi- 
Dirac distribution and the T^/n are the Unruh temper¬ 
atures. Now we can use 
calculate 


= (1+tan2(r/))iV^°2 to 

^BF 




(C15) 


"♦s, = E Nt = (1+‘“"<v)) E 

n n 

where from now on we identify and N^± , and 

is the negativity calculated from (jC12l) . Again, 

^BF 

(n) 

one can see that each NW is bounded from above by 

^BF 

I 


where 



(C16) 


The further N^2 for n ^ 0 can be obtained analytically 
and the negativity N*± is obtained to be 

^ R IT 




,± 

BF 


= ^i± 

^BF 


4- -T-(n -I- sinh^(r) cos^(r/) (tanh^(r-) -|- tan^(r/)) -|- 

^ 4cosh^(r) ^ ^ ^ ^ 

n—1 ^ > 


— \J'n? + 2sinh^(r) ((n + 2) tanh^(r) cos^(r/) -|- nsin^(r-/)) + sinh^(r) (sin^(r/) — tanh^(r) cos^(r/))^). (C17) 


To obtain a condition for vanishing negativity, we have 
a look at (jC16() and realize that the condition for entan¬ 
glement can be written as 


n‘gn9 < 1- 


(C18) 


Finally, we calculate the negativity of state 4'^^,, 
where we again, for the sake of concreteness, consider 
The relevant part of the reduced partially trans¬ 


posed density matrix is of the form 


( cLT ° + V 

\ 0 tan^(rj)c^_f;^ / 


(C19) 


where 


1 


Ce = - cos^(r/) 


tanh^’"(r) 


2 ^ cosh^(r) 

ecos(ry).^^ 


n+1 


cosh^(r) 

2/„ \ 


TS^sin (ry)+tanh (r) 


(C20) 


Due to the structure of (IC19I) the negativity again has 
the form (IC13F Then fV^± is calculated to be 






























21 


oo 1 2n-2/ \ 

+ 1) sin^(^/) + sinh^(r)) +ncos^{rf)+ 

BF j 4 cosh r) 

n—1 ^ > 

—2(— cos^(r/) + tanh'^(r)( ^ ^ sinh^(r) sin^(r/) + ^(^ + 1)^ sin"'’(r/) + - sinh^(r))+ 

+ tanh^(r) cos^(r/)((^ + 1) sinh^(r) — ^ ^ ^ nsin^(r/)))^), (C21) 


where 


_ 

'I'bf 2 (Zg)^ 


To extend the considerations of Secs. imiiiYi and El 
in 2d Rindler space to this spacetime, we consider the 
(C22) wave equation for a massless scalar field ip that is given 
by n\ip = 0. In the near horizon limit (ID3I) we can write 


again gives an upper bound on all the and 

^BF _ ^BF 

an lower bound on N^± ; see Fig. I8dl 

^ n vr 


Appendix D: Near horizon limit for a Schwarzschild 
black hole 


{Or - As2)ip{r],p,(p,e) =0, (D7) 

where (p, 0 are angular coordinates, is the 

d’Alembertian of 2d Rindler space, and A52 is the Lapla- 
cian of the 2-sphere. We are looking for solutions of the 
form 


In the presence of a Schwarzschild black hole the 
spacetime outside the black hole is characterized by the 
Schwarzschild metric, 

ds^ = ^1 — dt^ — ~ ^ (Dl) 


'(piV, P, = AadiV, p)lpang{(p, d) , (D8) 

that satisfy 

ORipradiil, p) =0, (D9) 

As^1paug{<p,e) =Q. (DIO) 


where G is the gravitational constant, M is the mass of 
the black hole, Rs = 2GM is the Schwarzschild radius 
and dil^ is the line element of the unit 2-sphere. In order 
to obtain the limiting form of m close to the horizon 
of a Schwarzschild black hole, we consider an observer 
placed at r = To with proper time 77 = (1 — 2GM/ro) t. 
Introducing 


The solutions of (ID9]) are the well known solutions of 
the Klein Gordon equation in Rindler space that we used 
above. The eigenfunctions of A52 are given by the spher¬ 
ical harmonics 0). The eigenvalues are l{l + l). So 

we pick the eigenfunctions with I = 0, i.e., 

Ipangi^, 0) = (cos(0)) = 1, (Dll) 


p^ = 8GM (r - 2GM ), (D2) 

we obtain the following metric up to terms of order 
in the near horizon limit 


ds — H — ^^2’ 


(D3) 


where 


16G2M2 

dsl ={2GMfd9?. 


2GM 

ro 


-1 




dp^, 


(D4) 

(D5) 


So (IDll) reduces to the product of two-dimensional 
Rindler space (ds|j) and a 2-sphere of radius 2GM {ds\). 
Comparing to the two-dimensional Rindler metric 
(ED, we see that the acceleration a experienced by an ob¬ 
server at fixed position r = ro is given by 


1 

4GM 



2GM 

ro 



(D6) 


where the P™(cos(0)) are the associated Legendre poly¬ 
nomials. We conclude that for the choice I = 0, i.e., 
zero angular momentum, we can describe the near hori¬ 
zon limit by restricting our considerations to 2d Rindler 
space. 

Therefore, in the following we restrict ourselves to wave 
functions ip of vanishing angular momentum satisfying 
(ID8II . When we consider maximally entangled fermion 
states dm) to hover over a black hole at some distance 
d = tq — Rs from the horizon, the system can be de¬ 
scribed in 2d Rindler space (for some more details on 
this correspondence, see M)- The analog of the Rindler 
vacuum |0)/ is the Boulware vacuum |0) b and the Unruh 
vacuum |0)f/ corresponds to the Hartle-Hawking vacuum 
\0)h. Further, the physical effect that causes the degra¬ 
dation of entanglement is now the Hawking effect. Near 
a black hole of mass M the acceleration a in (ED is set 
by dml). So, we see that the limit of infinite acceleration 
corresponds to the limit of tq approaching Rr. Consider¬ 
ing an observer stationary at a radial distance of rp, one 
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can write the acceleration parameter as 


r =arctanh(e ^ V 

(D12) 

rf =arctan(e ^ V '•o ), 

(D13) 


where ujg and fig are the rescaled frequencies ujg = 


AttRsuj and fig = AirRs^t. Plugging the acceleration pa¬ 
rameters (ID12I1 and (ID13p into the expressions for the 
negativities we obtained above, one obtains the negativ¬ 
ities of the respective states in the case that the acceler¬ 
ation is due to the presence of a black hole. 
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